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1 Introduction 



Effective time evolution of a quantum system is usually described by a dy- 
namical semigroup: a semigroup of completely positive, unit preserving, 
transformations acting on the algebra of observables of the system. A general 
form of generator of a norm-continuous semigroup was published in 1976 in- 
dependently by Gorini, Kossakowski and Sudarshan (for matrix algebras) 
on one hand, and by Lindblad 0] (for more general, norm-continuous case) 
on the other. It is usually referred to as the Lindblad form; it reads: 

A = z[H, A]+Y: V:AV^ - \{K A} (1) 
where H = H* is the Hamiltonian, { , } stands for anticommutator, and 

A = E (2) 

a 

In a contrast to a pure unitary evolution that describes closed systems 
and which is time-reversible, the second, dissipative part of the generator 
makes the evolution of an open system irreversible. This irreversibility is not 
evident from the very form of the equation, it is connected with the posi- 
tivity property of the evolution. Formally we can often solve the evolution 
equation backward in time, but positivity of the reversed evolution will be 
lost. 

We can also look at the dual time evolution of states rather than of observ- 
ables. For states, described by density matrices, we get: 

p = -t[H, a] + y: v^pv: - Up, A}, (3) 

where the duality is defined by Tr {Ap) = Tr {Ap). 

Here again only propagation forward in time is possible, when we try to 
propagate backward, then we will have to deal with negative probabilities. 
This irreversibility is reflected in the fact that pure states evolve into mixed 
states. How do mixed states arise? In quantum theory, similarly as in the 
classical theory, they arise when we go from individual description to ensem- 
ble description, from maximal available information to partial information. 
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Or simply, they arise by mixing of pure states. Pure states are represented 
by one dimensional projection operators P. If dfi{P) is a probabilistic mea- 
sure on pure states, then the density matrix p defined by p = / Pdfi{P) 
is a mixed state, unless dfi{P) is a Dirac measure. Contrary to the clas- 
sical theory, however, in quantum theory decomposition of a mixed state 
into pure ones is non-unique. So, for instance, the identity operator can be 
decomposed into any complete orthonormal basis: / = I], |i >< i\, thus in 
indenumerably many ways. This mysterious and annoying non-uniqueness of 
decomposition into pure states in quantum theory can be simply taken as an 
unavoidable price for our progress from classical to quantum, as a fact of life. 
And so it was. Yet it started to cause problems in quantum measurement 
theory. 

The first attempt to give a precise mathematical formulation of quan- 
tum measurement theory must be ascribed to John von Neumann. In his 
monograph ||^ he introduced two kind of evolutions: a continuous, unitary 
evolution U of an 'unobserved' system, and discontinuous 'projections' that 
accompany 'observations' or 'measurements'. His projection postulate, later 
reformulated by Liiders for mixed states, is expressed as follows: 

'if we measure a property E of the quantum system, and if we do not 
make any filtering which depends on the result, then as the result of this 
measurement the system which was previously described by a density matrix 
p switches to the new state described by the density matrix EpE/Ti EpE.^ 

A whole generation of physicists was brainwashed by this apparently pre- 
cise formulation. Few dared to ask: who are 'we' in the phrase 'if we measure' 
[^], what is 'measurement' 0, at which particular instant of time the re- 
duction takes place? How long does it take 0, if ever [||, to reduce? Can 
it be observed? Can it be verified experimentally 0, [ll|? Nobody could 
satisfactorily answer these questions. And so it was taken for granted that 
quantum theory can not really be understood in physical terms, that it is a 
peculiar mixture of objective and subjective. That it is about 'observations', 
and so it makes little or no sense without 'observers', and without 'mind'. 
There were many that started to believe that it is the sign of new age and the 
sign of progress. Few opponents did not believe completeness of a physical 
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theory that could not even define what constitutes 'observation' but 
they could not change the overall feeling of satisfaction with successes of the 
quantum theory. 

This situation started to change rapidly when technological progress make 
it possible to make prolonged experiments with individual quantum systems. 
The standard 'interpretation' did not suffice. Experimenters were seeing with 
their own eyes not the 'averages' but individual sample histories. In particu- 
lar, experiments in Quantum Optics allowed one to almost 'see' the quantum 
jumps. In 1988 J.R. Cook [|l^ discussed photon counting statistics in fiuo- 
rescence experiments and revived the question 'what are quantum jumps?'. 
Another reason to pay more attention to the notion of quantum jumps came 
from the several groups of physicists working on effective numerical solutions 



of quantum optics master equations. The works of Carmichael ||T3[, Dal- 
ibard, Castin and M0lmer [0, |15[, Dum, Zoller and Ritsch |jl6[, Gardiner, 
Parkins and Zoller Jl^, developed the method of Quantum Trajectories, or 
Quantum Monte Carlo (QMC) algorithm for simulating solutions of master 
equations. It was soon realized (cf. e.g. |18|, [19|, ^ |22|) that the same 
master equations can be simulated either by Quantum Monte Carlo method 
based on quantum jumps, or by a continuous quantum state diffusion. Wise- 
man and Milburn p3| discussed the question of how different experimental 
detection schemes relate to continuous diffusions or to discontinuous jump 
simulations. The two approaches were recently put into comparison also by 
Garraway and Knight [Q. There are at present two schools of simulations. 
Gisin et al. tried to reconcile the two arguing that Hhe quantum jumps 
can be clearly seen' also in the quantum state diffusion plots. On the other 



hand already in 1986 Diosi proposed a pure state, piecewise determin- 
istic process that reproduces a given master equation. In spite of the title 
of his paper that suggests uniqueness of his scheme, his process although 
mathematically canonical for a given master equation - it is not unique. 
This problem of non-uniqueness is especially important in theories of gravity- 
induced spontaneous localization (see |^ , also |2^, ^ and references therein) 
and in the recent attempts to merge mind-brain science with quantum theory 



31| , |32|| , where quantum collaps plays an important role. 
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In the next section we shall see how the situation changes completely with 
the new approach to quantum measurement developed by Ph. Blanchard 
and one of us (see |^ and references there )3. In Sec. 2 we will sketch the 
main idea of the new approach. We will also indicate infinitesimal proof 
of uniqueness of the stochastic process that reproduces master equation for 
the total system, i.e. quantum system+classical apparatus. In Sec. 3 we 
give concrete examples of non-unicity when only a pure quantum system is 
involved - as it is typical in quantum optics. In Sec. 4. we will give a rigorous, 
global proof of unicity of the process, when classical apparatus is coupled in 
an appropriate way to the quantum system. Conclusions will be given in Sec. 
5. There we also comment upon the most natural question: we all know that 
every apparatus consists of atoms - then how can it be classical? 



2 The formalism 

Let us sketch the mathematical framework of the "event-enhanced quantum 



theory". Details can be found in |^3[. To describe events, one needs a clas- 
sical system C, then possible events are identified with changes of a (pure) 
state of C. One can think of events as 'clicks' of a particle counter, changes 
of the pointer position, changing readings on an apparatus LCD display. The 
concept of an event is of course an idealization - like all concepts in a physical 
theory. Let us consider the simplest situation corresponding to a finite set of 
possible events. The space of pure states of C, denoted by Sc, has m states, 
labeled by a = 1, . . . , m. Statistical states of C are probability measures on 
Sc - in our case just sequences > 0, Y.aVa = 1- 

The algebra of observables of C is the algebra Ac of complex functions on Sc 
- in our case just sequences /„, a = 1, . . . , m of complex numbers. 
We use Hilbert space language even for the description of the classical system. 
Thus we introduce an m-dimensional Hilbert space Tic with a fixed basis, and 
we realize Ac as the algebra of diagonal matrices F = diag(/i, . . . , /m). 

^ Complete, actual bibliography of the Quantum Future Project is always available 



under URL: http: / /www.ift.uni.wroc.pl^ ajad/qf-pub.htm 
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Statistical states of C arc then diagonal density matrices diag(pi, . . . ,Pm), 
and pure states of C are vectors of the fixed basis of Tic- 
Events are ordered pairs of pure states a ^ f3, a ^ f3. Each event can thus 
be represented by an m x m matrix with 1 at the {a, j3) entry, zero otherwise. 
There are w? — m possible events. 

We now come to the quantum system. 

Let Q be the quantum system whose bounded observables are from the al- 
gebra Aq of bounded operators on a Hilbert space Tiq. In this paper we will 
assume T-Lq to be finite dimensional . Pure states of Q are unit vectors in TY^; 
proportional vectors describe the same quantum state. Statistical states of 
Q are given by non-negative density matrices p, with Tr (p) = 1. 

Let us now consider the total system T — Q x C. For the algebra At of 
observables of T we take the tensor product of algebras of observables of Q 
and C: At — Aq<Si Ac- It acts on the tensor product Hq <Si He — ®^=i7ia, 
where Ha ~ Hq. Thus At can be thought of as algebra of diagonal m x m 
matrices A — (oa/j), whose entries are quantum operators: Uaa £ Aq, — 
for 13. 

Statistical states ol Q x C are given by m x m diagonal matrices p — 
diag(pi, . . . whose entries are positive operators on Hq, with the nor- 
malization Tr (p) = X^qI^ {pa) = 1- Duality between observables and states 
is provided by the expectation value < A >p= I]„ Tr (AaPa)- 

We will now generalize slightly our framework. Indeed, there is no need 
for the quantum Hilbert spaces Ha, corresponding to different states of the 
classical system, to coincide. We will allow them to be different in the rest 
of this paper. We denote Ua = dim{Ha). 

We consider now dynamics. It is normal in quantum theory that classical 
parameters enters quantum Hamiltonian. Thus we assume that quantum 
dynamics, when no information is transferred from Q to C, is described by 
Hamiltonians Ha ■ Ha — ^ Ha, that may depend on the actual state of 
C (as indicated by the index a). We will use matrix notation and write 
H = diag{Ha). Now take the classical system. It is discrete here. Thus it 
can not have continuous time dynamics of its own. 
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As in ||3^ the coupling of Q to C is specified by a matrix V = {gap), where 



g^p are hnear operators: gap : Ti-p — > Ti-a- We put g^a = 0. This condi- 
tion expresses the simple fact: we do not need dissipation without receiving 
information (i.e without an event). To transfer information from Q to C 
we need a non-Hamihonian term which provides a completely positive (CP) 



coupling. As in ||3^ we consider couplings for which the evolution equation 



for observables and for states is given by the Lindblad form: 



Aa = i[Ha, Aa] + Yl ala^pgpa - \{K, Aa], (4) 
p ^ 



or equivalently: 



-i[Ha,pa] +Y.3apppglp - ^{A„,p„}, (5) 

p ^ 



where 



!^a = Y.Qloc9Poc- (6) 

The above equations describe statistical behavior of ensembles. Individual 
sample histories are described by a Markov process with values in pure states 
of the total system. In |Q this process was argued to be infinitesimally 



unique. For the sake of completeness we repeat here the arguments. First, 
we use Eq. @ to compute Pa{dt) when the initial state Pa(0) is pure: 

Pa(0) = 5aojV'o >< Vol- (7) 

In the equations below we will discard terms that are higher than linear order 
in dt. For a = we obtain: 



Paoidt) = lipo >< ipo\ -i[Hao, \ipo >< ipo\] dt- 



(8) 



while for a ^ ao 

Pao (dt) = gaao | ^0 >< V^O 1 9aao dt (9) 

The term for a = ao can be written as 

Paoidt) = Paoli^ao >< ^"0!, (lO) 
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where 

^"0 = - — ^ ) .S" ^ ,.\ ',"„ ' (11) 

and 

Pao = 1 - A(V'o, ao)dt. (12) 
The term with a ao can be written as: 

Pa{dt) = Pa \lpa >< 1pa \ , (13) 

where 

Pa = \\gaao'4'o\\'^dt, (14) 

and 

\\9aaoWo\\ 

This representation is unique and it defines the infinitesimal version of a 
piecewise deterministic Markov process. 



3 Non-uniqueness in the pure quantum case 

In this section we will show on simple examples the nature of non-uniqueness 
in the pure quantum case. 

For simplicity let us consider a two state quantum system whose algebra 
of observables is equal to M2x2- Let Tt be a dynamical semigroup with a 
generator L given by 

L(p) = apa* — -{a*a , p}, 

where a G M2x2- 



3.1 Pure diffusion process 

First let us show that the time evolution determined by L can be described 
by a diffusion process with values in CP^ ||34|| . 
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Let a two component complex valued process ipt = ii't^^'t)^ (prime denotes 
the transposition) be given by the following stochastic differential equation: 

dipi = fi{i)t)dBt + gi{il)t)dt, i = l,2, 

where Bt is a one-dimensional real Brownian motion and 

aiii^t) = 5I(< ^* >* ~ - ^ < «* >t< a >t ■iPl 

M'^t) = - <a>t ipl 
j 

<V't|a|^t> ilJt\a*\ipt > 

<a>t= -r- , <a >t= -r- 

< Ali^t > < Ali^t > 

Moreover let us choose an initial condition ipo = {zq,Zq)/ such that |2;qP + 
I^Qp = 1. Because and Qi are continuously differentiable (in the real sense) 
on \ {0} so there exists a local solution with a random explosion time T 
(see for example [^). But 

dl^f = ^Pld^Pl + i^id^i + dii^l^Pl],, 

where [ipljtpllt is the quadratic covariation of tpl and ipl- Thus 

dm,^i]t = \M^t)\'dt, 

and so 

dim' = Ei^id^i + ^id^i) + umi'dt = 0. 

i 

It implies that T = oo with probability one and so our process is a diffusion 
with values in a sphere S^. Let us define a process Pt with values in one 
dimensional projectors by 

Pt = \iJt >< iJt\ = Y.ijiij{e,j, 

id 

where e^j form the standard basis in M2x2- Then, using the equation 

dirM) = (M + M)dBt + {m + g^^i + hh)dt 
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we obtain that 

dPt = [{a- < a >t)Pt + Pt{a*- < a* >t)]dBt - ^{a*a, Pt}dt + aPta*dt. 
Since Bf is a martingale then after taking the average we get 

dE[Pt] = aE[Pt]a*dt - ^{a*a, E[Pt]}dt 
Let us define a density matrix pt = E[Pt]. Then 

pt = apta* - ^{a*a, pt}, 

and so the average of the diffusion gives the quantum dynamical evolution. 
Finally, we show that pt = J Pit, xq, dy)Py, where Py = \y >< y\, xq = 
>< ipol and P{t, X, dy) is the transition probability of the described 
diffusion. By the definition, P{t, xq, T) is the distribution of the random 
variable Pt° such that Pq° = xq. It implies that for every bounded and 
measurable function / defined on CP^ we have 

E[fiPr)] = J fiy)Pit, xo, dy). 
Let us consider a function given by f{y) = Tr{APy), where A G M2x2- Then 
J Tr{APy)P{t, xo, dy) = E[Tr{AP,^<^)] = Tr{Apt). 

SoTr{A J P{t, Xq, dy)Py)) = Tr{Apt) for every A and thus pt = J P(t, xq, dy)Pi 
with po = ^o- 



3.2 Piecewise deterministic solution 

On the other hand it is possible to associate with the same quantum dynamics 
a piecewise deterministic process, as in the method of quantum trajectories 
Now the situation is more complicated, because, in general, we can not 
replace the Brownian motion by the Poisson process. We have to solve a 
stochastic differential equation for an unknown process {Nt, ipt)- 

diPi = fii^lJt-)dNt + gi{'ipt)dt, 
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where fi and gi are prescribed functions, together with the following con- 
strain: Nf is a semimartingale such that 

a) [N, N]t = Nt, No = 0, E[Nt] < oo for all t > 0, 

b) for a given nonnegative function A : ^ R the process 
Mt := Nt — lo ^{'ips)ds is a martingale. 

It is clear that Mt will be a purely discontinuous martingale. A con- 
tinuous, increasing and with paths of finite variation on compacts process: 
Jo ^ii't)ds is called the compensator of Nt. In our case due to assumption a) it 
is also the conditional quadratic variation of A^^ [Q. The functional X{ipt) is 



called the stochastic intensity and plays the role of the intensity of jumps. Let 
us recall that for the (homogeneous) Poisson process Nt — /q Xds = Nt — Xt is 
a martingale. From the assumption a) above we obtain that Nt is quadratic 
pure jump, its continuous part is equal zero and ANg = {ANgY, where 
ANg = Ng — Ng- so it is a point process. Let us emphasize that in general 
it is not an inhomogeneous Poisson process since its compensator would be 
a deterministic function equal to E[Nt] [^. So it will be the case only when 
the stochastic intensity is a deterministic function depending on t. 
Moreover [A^, t]t = as A^ is of finite variation on compacts. It implies the 
following symbolic rules 

(rfiV)2 = dN, dNdt = dtdN = 0. 

From assumption b) we get dMt = dNt — X{tpt)dt. Let J-'t be a a-algebra of 
all events up to time t. Because Mt is a martingale, so E[dMt\J-'t] = what 
implies: 

E[dNt\J^t] = X{tlJt)dt 

see g. 

Till now the operator a G M2X2 was arbitrary. A particular simple case 
is if we take 

a* = a = 
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. Then L{pt) = apta — pt and so the intensity 

< -0*1% > 

what imphes that Nt — ^t- Because there is no deterministic evolution (we 
do not have the Hamiltonian part and the jump rate is constant) so in this 
case we can put = ^2 = and /i(V't) = "0* - ^l, f2{i^t) = -0* - "0* 
as the probabihty of a particular jump depends on the difference between 
and ■01 . Thus we arrive at 

Using the identity d[ijj\ -0^]^ = fifjdNt we get that rf||'0|p = and dPt — 
{aPttt — Pt)dNt. Taking the average we obtain pt — apta — pt, since Nt — Xt is 
a martingale. The above stochastic differential equation admits the following 
solution 

_ 1 + (-!)"• ^. 1 -(-!)"■ 



It implies that 



1 


2 


1 - 




1 


2 


1 + 




2 



.2 i + (-i; 
'0 7, 



- . V l-(-l)^* 

Pt = xo + yo- 



2 

where xq = \ipo >< ipo\ and 7/0 = |0o >< 0o|, 0o = (a + a*)^/;o = {zq, Zq)/. 
If we take 

as it is usual in quantum optics problems, then we have 



HA) 



So we need a point process whose rate function is random and the situation is 
slightly more complicated. We have to use the more general method described 
at the beginning of this paragraph. 
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Let us start with calculating functions gi, which are responsible for the 
deterministic flow. They are obtained by taking the derivative of 

exp(—^sa*a)ipt 

II,,. n 



with respect to s and in the instant s = 0. So we get 

aii't) = ^{-a*a + < a*a >t)ipt- 

It can be checked that the only functions fi which lead to the Lindblad 
equation are of the following type: 

where : — ^ R is an arbitrary Lipschitz function. Let us point out that 
if we put e*^ = ipf then we can write / in a compact form 



see 



38[] , but it needs a careful interpretation because zero can appear in the 



denominator. Again by simple calculations we get that = and 

1 / 2\^i\'\^p^\' rMli^^l'-l^ir' 



dt. 



But dNt = dMt + \{ipt)dt so after averaging we get the quantum evolution 
equation for pt = E[Pt]. 



4 Global existence and uniqueness 

After analyzing a typical example of non uniqueness in the pure quantum 
case, here we will return to the general scheme as described in Section 2. 
Let Tt be a norm-continuous dynamical semigroup on states of the total 
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algebra At corresponding to eq. (^. We extend Tt by linearity to the whole 
predual space At*, which is equal to At, because the total algebra is finite 
dimensional. Let E denote a space of all one-dimensional projectors in At- 
Because At = ®a=r ^i^a x no) we obtain that E = UaCPa and so is a 
disjoint sum of compact differentiable manifolds (complex projective spaces 
in TCa)- We would like to associate with Tt a homogeneous Markov - Feller 
process with values in E such that for every x E E 

Z{P,)= [ P{t,x,dy)Py, (16) 

where P{t, x, dy) is the transition probability function for the process and 
?/ — > Pj, is a map which assigns to every point y E E a one-dimensional 
projector Py. This leads us to the following definition. 

Let Ai{E) denote a Banach space of all complex, finite, Borel measures on 
E. We say that a positive and contractive semigroup Ut : A4{E) Ai{E) 
with a Feller transition function P(t,x,T) is associated with Tt iff Eq.|l6| is 
satisfied. 

Let us describe this notion more precisely. Let vr be a map between two 
Banach spaces M.{E) and At* given by 

7r(/i) = / jj,{dx) P^ 

JE 

It is clear that tt is linear, surjective, preserves positive cones and ||7r|| = 1. 

Proposition 1. Ut is associated with Tt iff kern is Ut - invariant and 
Ut = Tt, where Ut is the quotient group of Ut by kerir. 

Proof. Let Ut be associated with Tt. It implies that 

/ P{t,x,dy)Py =Tt{P.) 

J E 

thus for any /io G fcervr we have 

/ {Utfio){dx) P^ = / / P{t,y,dx) fio{dy) P^ = 
JE JeJe 

[ Tt{Py)^io{dy) = Tt[[ fio{dy)Py] = 

J E J E 
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and so Ut^o G kern. Moreover Wfi E A4{E) 

Utirii^) = n{Uti^) = I {Utfi){dy)Py = 
Je 

/ / P{t, X, dy)ii{dx) Py = Tt[ ji^dx) PJ = Tt-n^n). 

JE JE JE 

Now let us assume that Ut = Tf i.e. V/i G A4{E) we have Ut-n^ii) = Tt7r(yu). 
Let us take fi = S^- Then 

f/i7r(4) = TT{Ut6,) = f {Ut6,){dy)Py = 

J E 

f f P{t,z,dy)S,{dz)Py= f P{t,x,dy)Py 

JE JE JE 

and Ttnifi) = Tt{P,) so Tt(P,.) = Jj, P{t,x,dy)Py. □ 

It means that to find Ut is to extend the semigroup Tt from Ai{E)/ker ir 
to Ai{E) in an invariant way. It should be emphasized that, in general, such 
an 'extension' may not exist or, if it exists, need not be unique. We show 
that in our case, under mild assumptions, the existence and the uniqueness 
can be proved. 

Let us write the evolution equation for states in the Lindblad form 

p = -i[H, p\+y: ycpVk - \{p, E y^yc), 
k ^ k 

where H = diag{Hi, . . . , Hm), = H* G Mijia x ria) and Vk satisfy the 
following assumptions: 

a) (Vfc)aa = for every k and a 

b) if for some fc, /, a, [3 {Vkjap 7^ and {¥1)0,13 7^ then = /Q 

Let A be a densely defined linear operator on C{E) with D{A) = C^{E) 
given by 

= J2 Caix)f{Xa) - c{x) f (x) +v{x)f, 

In general we can allow for a weaker version: {Vk)a[3 and {Vi)a/3 ^ ^ 3c G 
C : (yk)ap — c{Vi)ap, but this simply reduces to b) above by substitution {Vk)ap ■= 
v/l + |cP(t4)„^ and {Vi)a,i3 = for fc ^ 
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where x G CP„(,, c«(x) = Tr {P^Wc^^cM^*^^), W^^a = Ek{Vk)aoa e 
rr — , G CPq, and x — f (x) is a vector field on E such that 

It may be easily checked that v{x) G T^CPa = T^E. Because 

''^ Tr (P. exp[-t Ea^ao W^«0«W^ao J) 

is an integral curve for f , so we have that f is a complete vector field. 

Theorem 2. A is a generator of a strongly continuous positive semigroup 
of contractions St on C{E). 

Proof. A = Ai + A2, where {Aif){x) = J2aj^ao Caix)S^^f - c{x)5^f and 
A2 = V. It is clear that Ai is a bounded and dissipative operator. It is 
also a dissipation i.e. Ai{p) > 2fAi{f) for f = f. Because A2 generates 
a flow on E given by f{x) fiotix)), where gt{x) is the integral curve 
of V starting at the point x, it follows that A = Ai + A2 is the generator 



of a strongly continuous semigroup of contractions (see for example |46|). 
Positivity follows from the Trotter product formula, since both Ai and A2 
generates positive semigroups. □ 

Let P(t,x,T) denote the transition function of St. 

Proposition 3. P{t,x,r) is a Feller transition function. 



Proof. By theorem 2.8 of ||39|, and by conservativeness of P(t,x, F), it is 
enough to show that for every x G P and for any / G C^{E) such that 
f{x) = 0, f{y) < Vy G P we have {Af){x) < 0. Because / has a maximum 
at X so {A2f){x) = 0. Moreover, as x G CPag for some ao and Ca(x) > 
Va 7^ ao we have 
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{AJ){x)= ^ c„(x)/(x,)<0 



□ 



Now prove that that our process reproduces Tf. 
Theorem 4. Let {Utn){T) := P{t,x,r)^{dx) for ^ e M{E). Then Ut is 
associated with Tj. 

Proof. At first we show that Wx E E 

L(P.) = [A(P)](x), (17) 

where L is the generator of Tf, A is the generator of St and P : x —>■ P^. 

Let a; G CPag. In H = ©™=i'^a let us choose any orthonormal basis 
{ea,iJtr=i\'..Zc.^ which e«,j^ G 7^^. Obviously, for any P^ G At* 
< ea,ijL{p)\ei3^i^ > = for a 7^ /3 and the same is true for [yl(P)](x). So it is 
enough to evaluate the (/?, i^, j/3)-th matrix elements of both sides of Eq. ( ]T7| ) : 



<e;3.,|[v4(P)](a;)|e^,,, >= ^ Tr (P.W,„,iy, 



Tr {P.W^^aW*) 



- ^ Tr iP.W^,c.W:j < e^,,^\P.\e^j^ > + 
<ep,^\{-t[H^„P,]~{P,, Y: W^,^W:^J+P.Tr (P, Wa,aW:j)\ep,^ >- 

a.^a.Q a.^a.Q 

^{P., E W^oaW:^a})\ep,,, > (18) 
On the other hand the /3-th component of L{Px) 

k k,a 

= W:^pP^W^,p + 5^,p{-i[H^,,P^] - i{P., Y Wo.,.W:j) (19) 
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Where the last equahty holds owing to assumptions a) and b) above. Taking 
the i/j, j^)-th matrix element of (|^) we see that it coincides with ([T8|), 
thus, due to arbitrariness of (/?, j/3), we have proved Eq. ([T7|). 
Let F denote the finite dimensional space of functions generated by a; — >< 
i)\PM >. It is clear that F = {f : f{x) = Tr{AP^), A e At}- So 
dimF = dim At- We show that F is the null space for kern. Let /(x) = 
< 'ipilPxli'j > and let /io G kern. Then 

Aio(/) = J fio{dx)f{x) =J2<^pi\ J fio{dx)P^\^pj > = 

Moreover because {A < ilJi\P\ipj >){x) = < ipi\L{Px)\ipj > we have that 
A : F ^ F and so St '. F —>■ F. It implies that Ut : ker n — > ker tt since 
Utfi{f) = fi{Stf). Let Ut be the quotient semigroup. Then 

limi[f/,(P.) - =limi[7r(f/A.) - = 

lim( J J P{t, z, dy)5^{dz)Py - P,) = (AP)(x), 

so Ut and Tt have the same generator and thus coincide. By Prop. 1 Ut is 
associated with Tf. □ 

We can pass to the uniqueness problem. Let us consider a class of 
Markov processes associated with a general nonsymmetric Dirichlet form S 
on L'^{E,dm) (here dm\cp^ is a positive f/(?^Q,)-invariant measure on Borel 
sets on CPa) given by the closure of: 

£{u,v) = / T{du,dv)dm+ / u{X.v)dm+ / (Y.u)v dm + / uvcdm+ 
JE Je Je Je 

+ f {u{x)-u{y)){v{x)-v{y))J{dx,dy) (20) 

for u,v e C^iE). In (H) we have : 

• X,Y - smooth vector fields on E 

• c e C°°(E) 
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• T - smooth (2,0)-tensor field, positively defined: T{du, du) > for any 

• A = {{x,x) E E X E} - the diagonal 

• J{dx, dy) - positive symmetric Radon measure on E x E\A satisfying: 

- lExE\Aiuix) -u{y)f J{dx,dy) < oo for any u G 

— the Radon derivative ■^'f^fy^ exists and is a Borel measure 

dm(x) 

• for any u > hold: /^(cm + X.u) dm, /_b(cu + Y.u) dm > 

It is worth to emphasize that such Dirichlet forms contain jumps, determinis- 
tic flows and diffusion processes as well. A straightforward calculation leads 
to the following result: 

Theorem 5. The generator B of the Dirichlet form S defined by ( pOD is 
given, in some coordinate system, by the formula 

{Bu){x) = U'u){x)+Y,V'{x){diu)ix)+Y,T^{x){d,dju){x)- f fi{x,dy)u{y) 

J E 

I ij 

(21) 

where / G C°°{E), V - smooth vector field on T - smooth, positive 
(2,0)-tensor field and fi - family of Borel signed measures, u G C°°{E) (the 
domain of B comes from closing {B,C°°{E)) ). The detailed form of /, V 
and fi is given by: 

/ = j^X.M + d,X' - c 
V = ^T^{d,M)d, + {d,T^)d, + X-Y 



and 



dm{x) dm{x) 
I 



dm{x) 

where we have used the following notation 
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• M is a coordinate of the volume form: dm{x) = M{x)dx^ A ... A dx^'^ 
for X e CP'' 

• 5x{-) - ei measure concentrated in {x} 

• an arrow indicates the variable the integral is evaluated over 

Remcirk 1. The generator B may be written in a fully invariant way: 

{Bu){x) = -cu ~ + ^ LiT(du)+uX] dm 
where Lzcu means the Lie derivative of the form cu, associated to the vector 



Remark 2. The proof of the Theorem 5 is straightforward - one should use 
the Stokes theorem for compact oriented manifold for the form 

oi = vi[T(du)+uX] dm 

and because E is without boundary {dE — 0) we have J^da — 0. Evaluating 
da we obtain the form of B. 

Using Theorem 5 and the property that Tr [L(Pa;)] = G -E we 
conclude that -B(l) = 0, 1-denotes the constant function taking value 1, and 
so B can be written in the following form: 



where {T'^^{x)) form a positive matrix and fiQ{x,dy) is a positive measure 
such that iio{x,{x}) — for every x & E. Its domain D{B) consists of 
C^-functions. 

Lemma 6. {Vk)aa = ^ Va G {1, . . . ,m} \/x,y G CPa such that Px-LPy 
the equality Tr[PyL{Px)] = is satisfied. 



field Z. 



{Bu){x) ^J2T''ix)ididju){x) +J2V'{x){diu){x) + 




(22) 
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Proof, let x,y e CPa and Px^Py Then 

Tr[PyL{Px)] = -iTr{Py[H^, p,]) + ^ rr[p,(y;p,yfe) J - 

k 

WTr[Py{Px, {VkVCU}] = J2Tr[Py{V:PxV,U] 
^ it k 

But 

SO the assertion follows. □ 



Wc arc now in position to show that the diffusion part is necessarily zero. 
Theorem 7. T^^{x) = for every i, j. 

Proof. Because 

B[Tr(PyP)]{x) = Tr[PyL{Px)] 

so, by the above lemma, for every a and every x,y E CPa such that Py^P^ 
we have that B[Tr{PyP)]{x) = 0. Let us denote the function z Tr{PyPz) 
by fy{z). Then 

{Bfy){x) = / ^Ji,{X,dz)^y{z) +Y^r\xMdjfy){x) + V^X^ fy) {x) 

ij i 
It is clear that fy is a smooth function and possesses a minimum at point x. 
So = and we arrive at 

lio{x,dz)fy{z) +^T^{xMdjfy){x) = 

But {didjfy{x)) and (T^^{x)) are positive matrices so, by Schur's lemma, 
{T^^{x)didjfy{x)) is also a positive matrix. It follows that 

Y^r^ix)d,djfyix) = 

ij 

Now let us introduce a chart at point x, let say, x — [(1,0,. ..,0)], (C/q, 0o) such 
that 

Uo = {[(^0, zi,..., Zn-i)] : ZiEC, \zi\'^ = 1, 7^ 0} 
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(f)o[{zo,Zi,...,Zn-l)] =( — ,...,^^) = {Xi,yi,...,Xn-l,yn-l), 

where Xi = -^e^, yi = Im^. Then 0o(a:^) = e R^^""-^). Let us choose 
= [(0,1,0,. ..,0)]. It is clear that Py-LPj. and so 



But for every j > 2 we have 



dx] 



(0) = 



hm i[^^^^(0, = 0, 0) - ^(^^(0)] = 

In the same way we prove that for every j >2 

^%^(0) = 
By positivity of the matrix D^{fy o (f)^^)(0) we obtain that 

Let A be an embedding A : CP^ CP^ given by 

A[(^o,^i)] = [(^0,^1,0, ...,0)] 

It is clear that x = A(n^) and y = A(n) for some unique tiq, n G CP^ = S'^. 
Let ipo be a chart at no given by 

■00 : CP^ - {n} C, ipo{rn) = p o 0o o A(rn), 

where p = C" — > C is the projection onto the first coordinate. So we may 
write that 

a K) (0) + 2a (no)^^^^(0) + 
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dqi 

where a^^(rfo) = T^^^{x), a^^(rfo) = T^lji^), a^^(n^)) = and qi{m) = 

Xi(A(m)), q2{'fn) = ?/i(A(m)). Let us change the chart ipQ onto spherical 
coordinates {9,ip), < 9 < n, < ip < 27i in such a way that 9{no) = n/2, 
'^{no) = i.e. no = (1,0,0) and ^(n) = 7r/2, (p{n) — tt i.e. n— (—1,0,0). 
Because 

/^j(m) = Tr{PnPrn) = <n,m>) = ^(1 - sin^cos(/?) 

so 

9^(^o) = 0, ^(no) = ^(no) = - 

which imphes that a}^{no) — a}'^{no) — 0^^(77^) = 0, where a'^ are the 
coefficients in the chart {9, (f). But it is equivalent to 

rii(^) = rii(x) = Tilix) = 

Changing 7/ = [(0, 1, 0, . . . , 0)] into y = [(0, 0, 1, 0, . . . , 0)] we obtain that 

and so on. Thus, by the positivity, T^^[x) = for every j, k. Because x was 
arbitrary the assertion follows. □ 

From the above theorem we conclude that the generator B has to be of 
the following form 

Bu{x) — V{x)u + / IjLo{x, dy)u{y) — Ho{x, E)u{x) 

J E 

with domain D{B) — C^{E) as is a closed operator. To proceed further 
we first need a lemma. 

Lemma 8. Let X be a tangent vector to CPq, at point P^- Then Px + X > 
O^X = 0. 

Proof. Because X e T^CPa so P^X + XP^ = X. It imphes that 



22 



PxXP^ = and P^XP^ = 0, where P^ = I - P^. It means that in a 

X* 
X 



X* 

basis PxH ® -P^^ form I I . So Pj; + X is a positive 



matrix if and only if X = 0. □ 
Theorem 9. B ^ A. 

Proof. Because A and B are generators of semigroups which are associated 
with Tt so for every x e E we have that [{B - A)P]{x) = 0. Let x e CPao- 
Then 

m . 

V{x)P fJ^oA^^dy)Py - iio{x,E)Px - 

E Caix)Px^ + c{x)Px - V{X)P = 0, 

where fiQ_a{x,dy) denotes the restriction of fiQ{x,dy) onto CP^. It is an 
operator valued equation so it has to be satisfied for every a separately. So 
for any a ^ ao we get 

/ IJ'0,ai^^dy)Py = Ca{x)Px„ 

JCPc 

which implies that iiQ^a{x,dy) — Ca{x)6{xa){dy). For ckq we have 

/ i^o,ao{x,dy)Py - i^o{x,E)Px + c{x)Px +V{x) - v{x) = 

Let US introduce a (x) = c{x) — fj.o{x,E) and w{x) = V{x) — v{x). Then 
taking the trace of the above equation we obtain a{x) < 0. Let us assume 
that a{x) < 0. It implies that 

If 1 

/ Ho,ao{x,dy)Py = P, - -——w{x) 

The left hand side of the above equation gives a positive operator and w{x) e 
TxCPao so, by Lemma 8, w{x) = 0. Thus we arrive at the contradiction 
because iio,ao{x, {x}) — 0. So a{x) — and we obtain that 



/ /^o,ao(a;, dy)Py + w{x) = 
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Evaluating the trace we get that /io,ao(^5 CPq,;,) = 0. Because it is a positive 
measure so it vanishes on every Borel subset of CP^o- So w{x) = too and 
hence A = B. O 

Thus we have the uniqueness. In the proof above we used repeatedly the 
fact that our Hilbert spaces were finite dimensional. In an infinite dimen- 
sional case the problem is much harder and we have no rigorous result. Our 
intuition is shaped here only by the infinitesimal argument of Section 3. 



5 Conclusions 

We have seen that the special class of couplings between a classical and a 
quantum system leads to a unique piecewise deterministic process on pure 
states of the total system that after averaging recovers the original master 
Liouville equation for statistical states. Irreversibility of the master equation 
describing time evolution of ensembles is reflected by going from potential 
to actual in the course of quantum jumps that accompany classical events. 
That is all fine but a natural question arises: what is classical? There are 
several options possible when answering this question. First of all the theory 
may be considered as phenomenological - then we choose as classical this 
part of the measurement apparatus (or observer) whose quantum nature is 
simply irrelevant for the given problem. Second, we may think of superselec- 
tion quantities ^ as truly classical variables. Some of them may play 



an important role in the dynamics of the measurement process - this remains 
for a while just a hypothesis. It is to be noticed that Jibu et. al (cf. [^], 
especially the last section 'Quantum Measurement by Quantum Brain'puts 
forward a similar hypothesis in relation to the possible role of microtubules 
in the quantum dynamics of consciousness. 

Finally, a careful reader certainly noticed that in the formalism of EEQT one 
never really needs C to be a classical system. It can be a quantum system as 
well. What is important it is that the the Liouville evolution preservers the 



diagonal of C. Thus the end product of the decoherence program [^, ^ 



can be directly fed into the EEQT event engine. The uniqueness result above 
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- will be immediately relevant also for this case. 
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